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ABSTRACT 


Seismic  arrays  are  multichannel  sensor  patterns 
Insnersed  in  a  multi -dimensional  signal-noise  field  and  the 
analytic  problem  is  hence  analogous  to  that  of  radar  antennas. 

The  subject  is  thus  opened  first  by  a  review  of  antenna  theory, 
considering  questions  of  aperture  width,  antenna  resolution, 
and  of  optimum  design  criteria,  and  secondly  by  a  review  of 
spectral  theci^',  including  special  examination  of  the  Ross 
"time  gates".  The  general  optimisation  problem  for  multichannel 
data  leads  to  large  systems  of  normal  equations  of  Toeplitz 
form  (as  presented  in  previous  reports)  which  require  recursion 
solution  techniques  to  be  computationally  feasible.  Such 
techniques  are  elaborated  here  in  terms  of  pjoiynojTda Is  orthogonal 
on  the  unit  circle.  The  specific  seismic  array  problem  is  then 
considered  in  terms  of  plrjtne-wave-front  signal  and  noise 
contributions  plus  incoherent  noise,  and  details  of  the 
"velocity  filtering"  method  are  presented.  All  practical  array 
filtering  rests  ultimately  on  empirical  measurements  of  signal 
and  noise  properties,  especially  of  spectral  behavior*  Spectral 
estimation  from  finite  array  measurements  is  the  final  question 
considered,  including  relations  between  continuous  and  discrete 
aperture  functions,  and  the  tabulation  of  aperture  functions 
with  their  windows , 
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SSISrttC  ARRAYS  FOR  THE  DETECTION  OF  MJCLEAR  EXPLOSIONS 

Snders  A .  Robinson 
May  196‘i 


;i,  OF  ANTENNA  THEORY 

Im.ys  of  detectors  have  velocity  discririlnation 
proper- ieCi  0.  hence  directional  sensitivity.  In  the  first 
■  t  iris  treatment  we  would  like  to  review  the  general 
theory  c*'  tne  directional  properties  of  antennas  In  order  to 
brlnt  ;nome  of  the  design  problems  for  arrays  of  specified 
direv  tl\  Ity  , 

An  antenna  may  be  viewed  as  a  spatial  filter,  and  so 
iv-s  a  bandwidth  that  is  determlnded  by  the  aperture  extent; 
therefore,  it  will  reproduce  only  a  finite  number  of  the  space 
harmonics  representing  a  desired  spatial  pattern.  From  this 
standpoint,  the  antenna  resolution  is  limited  by  the  highest 
space  harmonic  within  the  bandwidth  of  the  spatial  filter;  this 
bandwidth  in  turn  is  determined  by  the  aperture  size.  Never¬ 
theless,  this  well-known  limitation  on  antenna  resolution  may 
be  overcome  by  the  use  of  correlation  type  processing  of  the 
antenna  signals.  The  spatial-frequency  bandwidth,  and  hence 
the  angular  resolution,  of  an  antenna  system  thus  depends  not 
only  on  the  aperture  extent  but  also  on  the  time-frequency 
bandwidth  of  the  received  signal.  This  is  to  be  expected, 
since  the  aperture  extent  is  only  uniquely  defined  in  terms 
of  wavelengths  and  hence  the  signal  bandwidth  should  play  a 
part  in  determining  the  spatial  filter  characteristics. 

For  example,  let  us  look  at  a  radar  antenna.  Suppose 
that  a  linear  antenna  is  constructed  so  that  the  radiated 


electric  field  across  its  face  Is  sinusoidal  in  time  with  an 
amplitude  and  phase  depending  on  position  according  to  the 
complex-valued  function  of  position  A»C^)  .  The  superposition 
of  the  contributions  along  the  antenna  give  the  antenna  pat¬ 
tern.  By  examining  the  figure,  we  see  that  the  distance 
traveled  by  the  incremental  wave  from  the  position 
varies  with  the  position  X  . 


If  ^,CX)  ^  ,  then  the  sine  wave  r 

pos  ition  C  X ,  X  'Vdx^  is 

H  Co^tcOcCt-t-j+f 3^^ 

-  n  +  4^-^  j 


ecclved  froin  the 


where 


2TTC 


(  C  =  velocity  of  light). 


Let  us  use  a  refererce  phase  of  ~ 


2rr 


when  at  a  distar>^e  r„  . 


Then  the  received  signal  can  be  written  as 


M  Cc5l«.t+^^^'.r,-r)]dlx. 


We  will  now  make  the  approximation 


X  tavx0- 


which  is  good  for  small  ^  and  large  R  .  By  "far  field"  we 


mean 


o  ^ 

that  n  is  so  large  that  can  be  neglected;  the 


important  criterion  is  for 


<<  X 


8 


Vfhere  4  Is  the  aperture  'A'idth.  Also  we  shall  let  »  s 

r,-r 


HeiiCe  the  received  signal  may  be  written 
M  (ta5tcot+<^T^XB3dX. 


Adding  sine  waves  of  the  same  frequency,  the  amplitude  and 


phase  of  the  total  signal  -eceived  at  an  angle  %  is  given  Dy 


1 


Letting 


c0  = 


ITTX 

X 


and 


we  have 


aC8)  = 


I 


ni 

/  >■ 


IIT  / 

'isi 

X 


Acc^i) 


Thus  we  see  that  the  antenna  pattern  QLt8>  is  the  inverse 
Fourier  transform  of  the  Illumination  function  Aco5)  • 
Because  of  the  finite  aperture  width,  that  is, 

Acw5*0  {or 

the  antenna  pattern  is  the  Fourier  transform  of  a  band- 
limited  function.  Thus  a  narrow  beam  width  would  require 
a  wide  aperture  width  d  of  the  antenna. 


There  are  various  '^a.ys 


to  measure  the  beam  width  of 


an  antenna  pattcjrn,  \ 


oUJ^JX.'Ou  that  0-^^^  iS  r ea i  *  ihel 


1 

O'  - — — 


r 


( 

j  aifeidB- 


(where  §  ~  f  )  might  be  susll  not 

because  CXiBi  is  concentrated  around  on  the  C  -axis  but 

because  the  contributions  to  the  numerator  for  a>0  might  be 
cancelled  by  the  contributions  for  (k^C  .  Thus  we  see  that  it 
is  jd&\  that  is  involved  in  our  intuitive  notion  of  the 
spread  ofO-te)  on  the  &-axis.  Also  lCXie)(  allows  us  to 
consider  complex  as  well  as  real  0l(8^  .  Now  for 

analytic  reasons,  it  is  much  easier  to  work  with  instead 

of  10lC8>|  j  as  far  ar  spread  on  the  ^ -axis  is  concerned, 
is  as  satisfactory  as  fai8)|  >  although  obviously 
there  are  quantitative  differences  depending  upon  which  we 
use . 


Thus  the  beam  width  of  an  antenna  pattern  a(8)  may 
be  measured  by  the  quantity 


oO 


ot 


-■oO 


/OO 

&|alS)i  die- 

-»o 


10 


where 


' 

^  ^  AM 


It  is  always  worthwhile  to  consider  other  measures 
for  spread  along  the  ^-axis.  Another  measure  may  be  referred 
to  the  equivalent  rectangle;  specif ically,  the  measure  f  is 
the  width  of  a  rectangle  having  the  same  area  as 
having  the  same  peak  value  as  .  If  we  let 

denote  the  peak  value  of  >  then 


/  \m)y 


This  n^asure  is  not  good  for  antenna  pattems  for  which  fate) 
is  not  reasonably  block-shaped. 

Let  us  now  derive  expressions  for  d  and  f  in  terms 
of  the  illumination  function  .  Because  we  can  choose 

our  origin  of  coordinates  as  we  like,  we  may  assume  that  ^ 
and  9b  are  equal  to  zero.  This  involves  replacing  d^e)  by 
ate-6)  or  ate-e,)  as  the  case  may  be;  in  turn,  the  illumina¬ 
tion  function  is  modified  by  a  linear  phase  term  or 

respectively. 


Parseval’s  theorem  states  that 


oO 


arid  also  we  have 


a(e)^~j  j  Aci^)cico 

^  n  /  ^  J^ao 


.»oC 


Hence  the  measure  of  spread^  J  ,  is 


lI  ( 


hiO 

|Acd5)l^clc^ 


/  o<? 

r~  A(C1^ 


dul 


'-oO 


1  > 


the  width  of  the  equivalent  rectangle.  To  obtain  a  formula 
for  d}  in  terms  of  PioS)  ,  we  must  use  two  applications  of 
Parseval’s  theorem,  one  for  dUb)  and  the  other  for  ^0(0)  . 
Because 


ac0) 


-vaC 


we  have 

oO 

8  a(0)  dB- . 

qC 


That  is 


A  ftcu)  B  0.(6) 


is  a  Fourier  transform  pair.  Thus,  using  Parseval*s  theorer 
for  this  pair,  we  have 


,  oO 


1  dou. 


Thus  the  measure  of  spread. 


is 


oc 

y 

f  I 


These  formulas  for  f  and  allow  us  to  study  the  dependence 
of  beam  width  upon  fearturcs  of  the  illumination  function 


The  illumination  function  may  be  written 


Acus)  = 


where  Mc&i)  is  the  ma^iitude  and  (J’toi)  is  the  phase.  We  now 
want  to  show  that,  for  any  fixed  magnitude  function  Htos)  , 
the  phase  function  that  minimizes  the  beam  width 

is  a  constant  (or  linear  function  of  Od  ) .  Since,  under 
the  assumed  conditions. 


j  I  ^  ^ 


*  "-“OQ 


^oO 


1 


2. 


n 


IXed , 


t 


r-iS  ot* 
. *•••».  u-  M* 


Beca  us  e 


j  j  rie4diQ 

y~af} 


j  I  fi  6’'^ 


/ 

=  j  M  Jtb 


we  see  that  we  get  equality  only  if  4^  is  a  constant,  which 
proves  our  assertion.  Let  us  now  look  at  the  bean:  width 
Because 


=  (i  +  M(cd)  i 


we  liave 


3(^  »• 


y  frO 

J  Medico 


» *  I . 

because  the  factor  has  riO  effect  on  \  p‘  \ 

Because  H  and  are  real,  tfie  only  way  to  minimize 

|{q*4.^  for  fixed  M  is  to  make  <^'=0  , 

vihich  means  that  <!>  would  be  a  constant.  Nevertheless,  since 
we  only  need  to  make  M  4'^=^  0  >  the  points  at  which  Hcui)  =  0  need 
not  be  points  for  which  4>tuii  =  0  .  Thus,  suppose  exist 

except  possibly  at  isolated  points  where  Hc\.O3-0  and  at 


14 


such  points  let  hCoi'i  be  at  least  continuous .  Then  any 
that  is  constant  between  adjacent  zeros  of  H  minlmlzea 
As  usual,  any  linear  phase  can  be  subtracted  from  the 

minimising  without  changing  the  beam  width  because 

the  effect  is  only  to  translate  along  the  0^~axis. 

Suppose  that  the  antenna  pattern  0.16)  is  real. 

Because 


,uO 

Acw)=  I  ' 


oico)  e 


die 


) 


we  have 


oo  /  ^ 

ate)e^«sde  =  j  aiw 

=  Ac-oo)  , 

or 

Thus  the  phase  condition  for  minimum  beam  width,  namely 
constant,  requires  that  4^<.c^)*0  ,  and  so  that 

HCtO)  is  a  real  even  function.  In  turn,  the  antenna 
pattern  is  an  even  function. 

Let  us  now  look  at  the  effects  of  the  magnitude 
of  the  illumination  function  on  the  beam  width  of  the 
antenna  pattern  0(u))  .  Because  the  beam  width  can  be  made 
arbitrarily  small  if  the  bandwidth  of  is  made  large 

enough,  the  appropriate  problem  is  to  minimize  beam  width 
when  the  aperture  bandwidth  is  limited.  Thus  we  need  a 
measure  of  the  spread  of  the  Illumination  function  Acui")  on 


15 


the  cx)  -axis  to  provide  a  numerical  measure  of  bandwidth, 
such  notion  of  bandwidth  is  the  radius  of  gyration 


A 

j  |Acu))i^^u) 


where  we  assume  that 


/  CjO 

Q  =  j  CO  I  =  0. 


This  last  restriction  only  amounts  to  a  translation  of 
by  an  amount  53  ,  corresponding  to  multiplying  CXcO)  by  , 

which  has  no  effect  on  Ct  (or  f  ).  Another  notion  of  band¬ 
width  is  the  natural  one  for  antennas  with  finite  aperture 
width;  here  it  is  assumed  that  outside  of  an 

interval  (corresponding  to  the  aperture)  on  the  cO-axis,  We 
define  “fe.  as  half  the  length  of  the  interval;  in  terms  of  our 
previous  notation 

A(dj)=;0  =  ^ 

where  (j[,  =  aperture  width  and  ^  (where  C  = 

velocity  of  light). 

We  now  wish  to  shov/  that  broad  illumination  band¬ 
width  gives  a  narrow  beam  width  (that  is,  a  good  antenna 
resolution) , 


ana 


,  We  have 


Prom  the  derinitlc=ns  of  (k 


J  I  lAi^Ju) 


By  use  of  the  Schwarz  inequality,  it  follows  that 

/  ^  2.  f  ^  f  ^ 

1/  QAA'dtoj  ^  J  1<^'AiMcoJ  lA^j^dlu) 

'  Aao  / 


SO 


oo 


J  <^AA  dcoj 

fo^ 


^  ot 


We  assume  that  the  phase  4^(fcd)K0  in 


/\(to)fc  (^cq) 


so  that 


Acw)  =  M(«)’) 


Is  real.  H  ^nce  we  see  that 


'  =J.  dft' 

2  dco 


thereby  giving 


for  the  numerator  of  the  bove  expressiun 


If  we  Inteerat! 


)y  parts,  we  therefore  obtain 


i  “*  J- 


H^dluO 


^2 

for  the  numerator,  and  since  cOA 
expression  reduces  to 


as  leal  ,  this 


Hence 


i  ^ 


Which  is  the  fundamental  expression  relating  the  measure  Ui 
of  the  beam  width  (or  antenna  resolution)  and  the  measure 
of  the  Illumination  bandwidth  (or  aperture  width).  Thus 
if  we  wish  to  narrow  the  beam  width  of  the  antenna  pattern 
it  is  necessary  to  broaden  the  band  width  of  the  antenna 
illumination . 

The  lower  bound  of  -j  for  is  actually  obtained 

for  the  Gaussian-shaped  illumination  function 


Acoii  =  C! 


(where  C  =  a  positive 
constant ) 


which  yields  the  antenna  pattern 


(which  is  also  Gauss Ian- shaped) 


y 

ace)^  C 


(where  C  s  a  positive 
constant) . 


Nevertheless,  there  is  no  upper  bound  for  . 

A  similar  result  to  can  be  obtained  by 

using  the  width  of  the  equivalent  rectangular  f*  to  measure 
beam  width  of  the  antenna  pattern  CXC^)  and  using  finite 
aperture  width  to  measure  the  bandwidth  -k.  of  the  illumination 
A'v  .  we  suppose  that 


A(Q)  =  0  for 


Hence  the  width  of  the  equivalent  rectangle  is 


The  integrand  of  the  denominator  is  A  ,  which  we  may  write 
as  1-A  in  order  to  apply  the  Schwarz  inequality.  Thus 
we  have 


-k  --k 


satisf led: 


Therefore  the  following  inequality  is 

if  >'TT. 


The  equality  is  satisfied  if 
is  a  constant  for  {c0|<  ^  , 
is 


and  only  if  the  illUjnination 
in  v/hich  case  the  antenna  pattern 


a(6)  =  C  - 


(where 


a  constant) 


In  view  that 


ty  -tt 

for  a  constant  illuniinatiori  function  over  the  finite  aperture, 
that  is,  for 


constant 

0 


for 

-or 


It  follows  that  this  Illumination  is  optimum  in  the  sense  of 
giving  the  smallest  possible  equivalent-rectangle  width  for  a 
given  aperture  size.  Nevertheless,  if  we  measui^e  the  beam 
width  (l.e.,  antenna  resolution)  in  terms  of  radius  of  g^^ration 
Instead  of  in  terms  of  equivalent  rectangle,  we  obtain  a  dif¬ 
ferent  result  for  the  optimum  illumination  function  for  a 
f-’nite  aperture. 


The  finite  aperturf^  4- 


Lun  function  satisfi 


,riction 
the  condition 


r*-t5  7-.  '+ 


he  ill 


AcoC))^0 


i  01 


l|  >  t 


The  result  that  we  seek  is  the  siiape  of 


^i^) 


for 


such  that  the  product 


Is  a  minimum.  To  minimise  d  ,  A  will  have  a  constant 
phase  that  we  can  take  to  be  zero,  so 

is  real.  We  recall  that 


Thus  we  wish  to  minimize 


F<X> 


(AVj.co 


■<30 


( A')^  doO 


>  1  O  t  A 


A  dw  ^ 


Using  the  calculus  of 


>e 


vdi  4.0- 


ions ,  the  solution  comes  out 


s 


e  I  31  \ 


(where  (^  =  a  constant) 


cnc  the  minimum  of  Ck  is 


z 

d 


i 


^ie 


JZ 

f  a 


<.3t 

4 


4fi^ 


Thus  beam-width  parameter  ci^  and  aperture-width  parameter 

a  - 

—  satisfy  the  Ineaualitv 


A 


JL 

2 


with  equality  if  and  onl’ 

'1c  1  •  f'!  £5  •■4  O  »  » /a  W*.  ^  J  T  T  -  , ^  . 

4.ij  ouCfi  Oil  X j. xliiaxfiii 

antfsnrias  . 


if  a  cosine  illuminatj 
tion  function  is  in 


i  iUiiCua-OfI 


c ornmori  liso  ^ri 


hat 


iuniTiary^  we 


ooserve 


t  Vl 


■*.  ^  ^  W'  ks 


(1)  to  have  a  narrow  beam  antenna  i^attern 
It  is  necessary  but  not  sufficient  to 
have  a  broad  illumlnat ioii, 

(2)  phase  of  the  illumination  function, 

other  than  a  constant  or  a  linear  phase 
(which  shifts  the  antenna  pattern  with¬ 
out  ciianging  its  shape)  increases  the 
beam  width  of  the  antenna  pattern,  and 
there  is  no  upper  bound  for  this  increase. 


*0 


innij 


*  rr:*  ^  * 

rLCr  1. 


We  now  want 


’evle'w  time-f requeriC2/  aspects 


of 


signals.  Nfeny 


O  1  : 


theory  may  be  traced  to  Fourier  analysis 
found  to  be  extremely  well-suited  to  the 
cal  statement  of  many  natural  phenomena, 
acterlstlc  of  the  Fourier  method  is  that 
equivalent  statements  in  either  the  time 
quency  domain. 


which  has  been 
precise  mathematl- 
An  important  ctiar- 
it  gives  completely 
domain  or  the  fre- 


There  are  three  main  mo  .sis  for  the  Fourier  approac.h 
th.at  can  be  set  up,  namely 


{!)  the  Fourier  series  approach,  with  the 
attendant  assumption  that  the  time 
function  is  a  periodic  function  of 
continuous  time,  with  the  result  that 
the  spectrum  is  made  up  of  lines  a., 
discrete  frequencies  that  are  integer 
multiples  of  the  fundamental  frequency 
(equal  to  the  reciprocal  of  the  period). 

{2}  the  Fourier  Integral  approach,  with  the 
attendant  assumption  that  the  time  func¬ 
tion  is  an  aperiodic  integrable  fujiCtlon 
of  continuous  time,  with  the  result  tint 
the  spectrum  is  also  an  aperiodic  inte- 
gi’able  function  of  continuous  frequency. 

(3)  the  z-transform  approach,  with  the  atten¬ 
dant  assumption  that  the  time  function 
is  a  functlcn  of  equally-spaced  dis¬ 
crete  time,  with  the  result  that  the 
spectrum  is  a  periodic  function  of  fre¬ 
quency,  Hence  the  spectrum  is  a  Fourier 
series  with  the  time  function  as  its 
Fourier  coefficients. 


For  the  purposes  of  this  section, 
the  B'ourier  Integral  approach. 


model  (r): 


we  shall  use 
The  classical 


definition  of  the  Fourier  transforn:  is  as  follows:  Given 


.me~functlon  ,  then 


i  its  Fourier  transform  Is 


frequency-function  given  by 


,*o 


—  «c> 


where 


t,  =  time  (say,  in  seconds) 

O  =  radian  frequency  (say.  In  radians  per  second) 

It  is  necessary  to  state  various  further  remarks  about  -ftil 
and  FCcoi  ,  but  we  will  assume  that  the  reader  is  already 
familiar  with  them.  Making  use  of  the  Euler  identity 

-  COScOt-cSl^uSt 

we  see  that  the  Fourier  transform  is 


oO 

=  (cosine  transform)  -  -t,  (sine  transform). 

Much  of  our  discussion  will  be  concerned  with  the  cosine 
transform  for  a  parallel  discussion  would  apply  to  the  sine 
transform. 


Ft)  dlt 


It  is  often  the  case  that  only  a  finite  portion  of 
the  time  function  is  available.  As  a  result  an  error  would 
result  in  the  computation  of  the  cosine  transform  because  of 
the  use  of  a  finite  instead  of  an  infinite  interval  of  inte¬ 
gration.  This  error,  called  the  truncation  error,  can  be  of 
major  consequence,  and  so  v/e  will  discuss  various  auproaches 
to  this  problem  in  this  section.  The  effect  of  the  truncation 
error  on  transforms  Is  to  superimpose  a  relatively  large  am¬ 
plitude  ripple  upon  the  correct  transform;  this  ripple  is 
often  called  the  Gibbs  phenomenon.  Hence  we  will  look  for 
various  suitable  modified  computational  techniques  that  lessen 
this  spurious  ripple.  More  specifically  we  will  look  for  some 
suitable  function  that  we  sliall  call  the  time  gate 

associated  with  the  approximating  method. 

Denoting  the  Fourier  transform  operator  by  F  and 
the  approximation  operator  by  ,  we  have 

F(wd)^itn  = 

for  the  approximate  trarisform  F.cuj’i  . 

The  Fourier  transform  of  the  time  gate  is  called  the 

frequency  window  ;  that  is, 

V7c«^  =  F  C  i-ntd . 


Because 


Few-)  =  F 


we  can  find  the  approxlnkation  in  terins  of 

and  ,  We  use  the  familiar  rale  ti.it  the  transform  of  a 

prodact  is  eqaal  to  the  convolution  of  the  transforms  of 
each  factor  in  the  product.  Thus 

p(Xa)3 

A 

Where  %  denotes  convolution,  that  is 

/>o 

W(«o^*Fccd)=i  1  Vjcco-Gij'i 

Thus  we  see  that  the  approximation  is  obtained 

from  the  correct  transform  Fctb)  by  convolution  of  the 
correct  transform  with  the  frequency  windov^ ..  Let  as  now 
look  at  the  process  of  convolution  in  orde*  to  see  how  the 
approximation  differs  from  the  correct  transform.  We  see 
that  is  the  same  as  reversed  in  frequency 

direction  and  shifted  by  the  amount  CO  .  Taking  the 
product 

and  inte^ating  over  cOi  gives  the  value  of  the  approximate 
transform  at  the  frequency  cO  .  Hence  it  may  be  said  that 
we  "look  through"  the  frequency  window  (reversed  and  centered 
at  fid  )  at  the  correct  transform  in  order  to  obtain  the 
approximate  transform. 

The  frequency  window  is  Independent  of  the  time 
function  being  transformed;  Instead,  it  is  characteristic 
of  the  approximation  method  used  and  gives  a  complete  measure 
of  the  difference  between  the  correct  transform  and  the 
approximate  transform  for  any  time  function.  Ideally  we  would 


like  the  frequency  window  to  be  the  Dirac  delta  runction  for 
then  the  approximate  transform  would  be  identical  with  the 
correct  transform.  Practically  we  will  seek  a  frequency  window 

that  approximates  the  delta  function;  then  the  convolution 
process  for  a  given  frequency  averages  the  correct  transform 
over  a  narrow  band  of  frequencies  around  this  frequency,  with 
almost  no  contribution  from  frequencies  far  away  from  this 
given  frequency. 

^^et  ub  now  turn  our  attention  to  some  gates  and 
their  windows.  We  have  the  following  table: 


TIiME 

GATE 

FREQUENCY  WINDOW 

Box  car  of 

length  ZT  : 

Dirichlet  window: 

f  ' 

Vlt)  =  1 

■Jot  ItUT 

Vc«>  =  2t 

1  0 

H1>T 

=  irt) 

Triangle  of  length  2.T  * 


S/Sii)- 


1--^  -Joi-  Itl^T 

0  40’’  ltl>T 


Fejer  window; 

Z 


Delta  function  at  4,0- 

l  .  A  .X 


CoS  oist 


Von  Hann  gate: 

Wet)  * ! 

(  0  -fo*  \tl>T 


Richard  Hamming  gate: 


0.5H-40.4iC5»4^  ^ti^T 


V. 


0 


■(«t  \t\>J 


V3tt^  =  |.ct> 


=  CoS 


Ross  gate; 


I  0 


■^OT  1U<T 
4ov  ltl>T 


FREQUENCY  wIND<^ 


2  L 

Von  Hann  virndow: 


VTc^)  ^  Q^(  t- 


s:«ciT 


Richard  Hamming  window; 

*  0.l3  Q*C^>t  ^  )t  05^^ 
4oa3QoCc;>-^) 


WCo)= 


The  following  result  Is  known  from  the  tlieory  of  the 


Fourier  transform; 


If  a  time  function  has  ft-?  continuous 
derivatives  and  at  most  a  step  discon¬ 
tinuity  in  its  derivative,  then 

its  Fourier  transform  vfill  have  an 
envelope  of  order 


Let  us  use  this  result  to  compare  the  above  time  gates.  Hie 

th 

box  car  has  step  discontinuities  in  its  0-~-  derivative;  hence 
its  window  lias  an  envelope  of  the  order  of  *73*  ,  as  seen  by 
the  actual  expression  for  the  window.  The  trxangle  has  step 

g  fc 

discontinuities  in  its  1 —  derivative;  hence  its  window  lias  an 
envelope  of  the  order  ,  as  seen  by  the  actual  expression. 

The  Von  Hann  ^te  has  step  discontinu.Hles  in  its  second  der¬ 
ivative;  hence  its  window  has  an  envelope.  Hie  Ross 

gate  has  step  discontinuities  in  its  derivative; 

hence  its  window  has  an  envelope. 


It  is  interesting  to  consider  the  characteristics  of 
the  family  of  Ross  gates.  The  first  few  members  are  plot tec 
in  the  following  diagram: 


We  See  that  higher  values 
from  values  whose  abclssa 


of  i  means  less  and  leas  contribution 
Is  close  to  t’iT 


Let  us  now  find  an  interpretation  for  this  family. 
Let  us  consider  the  '"ase  where  the  function  to  be  transformed 
is 


|ct)  =  Cci 


for  some  arbitrary  frequency  cO  •  This  choice  has  the 
advan  age  that  we  know  that  the  normalized  correct  transform 
of  this  function  is  unity  for  co  and  zero  for  all  other 
frequencies.  The  problem  is  to  find  out  what  is  the  effect  of 
calculating  the  approximate  transform  given  by 


instead  of  the  correct  transform  given  by 


Um* 


T 

dlt 


) 


and  how  the  result  varies  for  T  and  c5»  .  The  frequency 

Q,  may  be  called  the  scanning  frequency , 


is 


For  the  case  when  ,  the  approximate  transform 


2.UJT 


2(*)T 


and  for  the  case  when  ,  the  approxinate  transform  is 


^  f  '  ^  dI,w.Ct3-tJ,)T  J  3-vv^C^=h 


•'  C‘ 


The  result 
figure. 


1  + 


luoT 


for  the  case  u}i  =  u5  is  plotted  in  the 


Figure  3.  Dependence  of  the  approxinatc  transform  on 
truncation  length  y. 


From  the  figure  we  see  that  the  approximation  oscillates  about 
the  correct  answer,  1,  and  gets  closer  and  closer  to  this 
correct  answer  as  T  is  increased.  The  problem  now  is  to 
find  a  way  to  utilize  this  oscillating  approximation  so  as  to 
find  a  more  refined  approximation.  It  appears  that  a  linear 
weighting  of  the  oscillating  estimate  shown  in  the  figure 
would  be  a  good  choice,  and  in  fact  corresponds  to  F'  jer 
summation  of  the  original  integral. 


We  thus  have 


I  -  tutioT 


) 


CoT)'-^ 


This  function  is  plotted  in  the  figure: 


Figure  4,  Dependence  of  the  second  approxlnate  transform 
on  truncation  length  T. 

For  most  T  values,  this  function  is  closer  to  the  correct 
valuer  1,  than  the  previous  approximation,  and  thus  the 
weighting  has  helped  matters  to  some  extent .  llie  envelope  of 
this  aporcximation  is  •  Nevertheless,  this  choice 

of  weighting  has  introduced  the  problem  that  the  approximation 
no  longer  osclllatec  about  the  correct  value 


Pur  the  third  approxlntition  let  us  choose 
weightlrig  function  31*/^’^  ,  thereb'/  obtaining* 


f 


3 


( 


lusT 


Cos  2olT  ] 


I 

This  function  not  only  has  the  envelope,  but  also 
the  approximation  still  oscillates  about  the  corx’ect  value, 
which  means  that  the  process  can  be  continued. 


For  the  scanning  frequency  tO.+uj  ,  this  approximation  is 


T" 


^  CojCO-cO 


.)t) 


C03u).)T 


which  has  the  same  order  envelope  around  the  correct  value  0, 


We  are  now  in  a  position  to  derive  the  Ross  class  of 
time  gates.  The  first  stage  gives  the  (unnormalized) 
approximation 


1 


The  second  stage  gives 


t, 


r  ^  r 

j  t.dt,  )  -f It)  Coi  oit  cit 


this  expression 


Interchanging  the  order  of  integration, 
becomes 

-fitKosCotdrt: 

so  the  second  stage  corresponds  to  the  weighting  function 
or  time  s?ite 

2 

Similarly  the  time  gates  for  higher  stages  can  be  found. 

It  is  instructive  to  look  upon  the  Ross  time  gates 
as  a  means  of  countering  the  effects  of  truncation 
bT  smoothing  the  approximation  at  the  price  of  losing  a 
certain  amount  of  frequency  information*  When  is  multi¬ 
plied  by  ,  the  effect  is  to  make  the  values  of 

for  large  tt|  not  as  important  in  the  approximation  as  those 
for  small  *t  .  Ihis  leads  to  smoothing  since  each  cuccesslve 

value  contributes  less  and  less  to  the  accunnjlated  Integral. 
Nevertheless,  if  we  take  a  time  gate  for  a  too  large 

value  of  4c  ,  then  we  essentially  reject  a  very  large  portion 
of  the  given  finite  record  of  -^tt)  ;  and  the  resulting 

approxination  achieved  by  using  a  more  moderate  value  of  -k 
Thus,  if  we  try  to  achieve  too  much  anxiothing,  we  may  sacrifice 
significant  frequency  information;  it  is  necessary  to  strike 
a  balance  between  these  two  conflicting  goals  Excellent 
empirical  and  theoretical  treatment  of  this  problem  is  given 
in  the  original  work  of  Ross  (195^), 


dt 
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■  i  uj  i 

*2  U  W  C 


re  wish  to  consider  properties  of  polynomials 
on  the  unit  circle  (Geronimus,  I960).  The  poly¬ 
nomials  are  defined  to  be  orthonormal  with  respect  to 

the  spectral  density  ^(j)  on  the  unit  circle  |2j  =  !  ;  that 

is. 


f  0 


V  I 


where  ~  VO., 

Also  let  us  define 


with  oi«>0 


4,h2)  =  2''fy^) 


so  <^*(o)=o(t 


Let  us  put 


fn-i 


“  +  . +  W„-^., 


. 


for 


and  then  multiply  each  side  by  4^(Z) 
set  ^  integrate. 


We  obtain 


,05 


y 


o.S 


,  o»?' 

'-o.S' 


since  <^,Ce'^’'‘b/s-^’'“"-'H  , 

,  ,  is  a  polynomial  in 

\  TL^T  integral  Is  zero,  and  by 

introducing  the  notation 


w;^  ^fr  i-0.1,2, 


■e  obtain 


and  hence 


•^*Z*'-2'->|’x*Vl^ 

Oln  -|{,0 


4> 

Thil 


Qln+t  -1=0 


_  2^ 

cU 


so 


By  Setting  i  — 0  wg  i  iriQ  »  Ihat  is^ 


d 


k*‘l 


SO 


Xh="'“ 


^■r\  1 1  o 
i.  2  i  L* 


>Vit! 


d, 


fHA-i 


'^i 


30  the  following  recursion  relation  holds; 


I 

In  this  relation,  let  us ’replace  2  by  ~  ,  thereby  obtaining 


= lilMl  ^ 

dft’t'i  div  d>\ 


and  multiplying  by  2 
obtain 


n4l 


and  taking  complex  conjugates  we 


^.A)  _  7"?,d)  ,  JJT)  z2•'+^i:) 


iK+i 


dn 


ln%l 
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which  is  the  recursion  relation 


and  define  the  parameters  Oh  by  the  forimila 

“T— - - 

dn^\ 

Then  the  recursion  relations  become 


Let  us  now  find  the  connection  between  the  parameter 
lAn  and  the  autocorrelation 

JL 

’’a 

The  roepiitz  determinant  is  defined  to  be 

^  t*  r,i 

^  =  lri.4(  »  '■»  r.  r.^,  Ah-I  • 

-  '  *  ) 


39 


and  it  is  seen  that  the  polynomials  can  be  written 


r,  r.,  r.2 


tin 


=  -aw. 


t  z  z* 


Let  ns  introduce  the  notation 


i 


where  and  la)  are  arbitrary  polynomials  in  Z  . 

Using  this  notation  we  see  that 


Ah- 1 


r.  r., 

V-,  r* 

^h-4  tii-2 
0.2“’  (2,2'') 


r-h 


n,  I 

CZ-Z")! 


I 


A,., 


K 

r., 

r.K 

r, 

ro 

1 

Ah 

*: - - 

t-.a 

Vmi 

t-l 

An-) 

^11 -t 

n> 

Prom  the  relations 


4u 


n 


it  is  possible  to  find  a  recursion  involving  only  ai 

follows.  Solving  for  we  find 

vt  ^ 

sr  -/fi- _ Zmi  . 

a." 


and  hence 


Which  gl'ves 


O.K 


tt4«2  =  c  i  -  ia„ih  tn**. 

Similarly,  we  find  that 

~  CXh4.i)  ^n+i^  (  |Qh|  ^ 

Prom  this  last  expression  we  obtain 


a.(t.  z**')  =  7"^') +aK«Ct.,  ,2"')  -  i 


which  is 


(l-lQ.l^)4r=0 

«n 


since 


Thus 


ii  IF  -J- 

iro  d 


.and 


c^n  dl  Cil  d?  ’dX 


This  last  step  follows  by  Induction  since 


-i—  =  /  -L  _ 

^  di."  d,^ 

_jQ/  la.i^ 

^ 


42 


lh  summary  we  have 


(1) 

(2) 

(3) 

(^) 


4>. 

C  ^  j  - 


Kronecker  delta 


f uncti 


(7) 


('-lQ.h(i-ia,|’). 


(I-IQh..^) 


di 


(Qtt-.{' 

oi«^ 


i 

(  n  r-i 

\ 

~  di»l  / 

r,  r* 

•  ^  « 

"for  ^1  =  0, 1,2^... 

\ 

r*  J 

> 

li 

7 

4 

! 

n. 

AW  1  A  ^  ^  " 


•f-  K  ^  Fns  1 1  ^  4  i-.  Q  r--r-  ^  ^  -»  —  o  -t.  A  ^  ^  1  j  _ _ _  ^  l _ 

viitj  i;i«Xi.--i-wii£inriei  cabe  iuxiuvsins  ^*1^ 


^  •■^  e^  .— »  ■€  •*  T  A  i  4-  "L—  _ 

i^iVeli  j.ri  tne 


It?  munugraph  Studies  In  Uptlmum  Filtering  of 


Sinsle  and  Wjltlple  StOf'Pa 


O  wOCrid^j  u  ic  Processes ,  We  now  use  matrix 


riutatlon,  we  nave  tPe  forward  polynonilal 


+  - i" 


C  i  ) 


and  the  backward  polynomial 


wnere 


-  bhn(x)  ’  '4  bhii  ?  )  +  b^o  C  b^o“  I  ) 


[  Q«o j‘” 


I  )r 


'  nvi 


6  y 


[  oi„,c,  ■■•,0,^;] 


with  d„  =  a„c  u., + ■  ■  • + a»n  r, 


and 


with  —  bnn^t  +  -  *  *  "I-  bno  . 

(Note,  this  is  a  new  use  of  the  symbol  ,  corresponding  to 
the  usage  in  the  above  monograph). 


where 

H-vl. 

Thu?  the  polynomials  satisfy  the  recursion  formulae: 


Ci 


rtVi 


n 

f  -■  i 


»V+i 


Ctn  + 


r,  A  ■'  i  A  Q 

Pr\=-i  I 


+  +iaipi  + 


i^n+t 


=;  3r\''K 


K 


These  ire  the  analogut^s  of  the  formulae  found  in 
the  single  channel  case.  Just  given. 


The  general  filter  polynomial  is 

^J.Z)  K  f  f  n,Z  4- 


t  In  n 

I  i  4U=i-  C 


r, 


H'Vi 


■*  ■w  h 


X^-\  -  .  fe 


h.,=  fn.''n-  +  -+t«'r, 


TV-.Jto 

i  i  i  ^  i5 


n  orrt.»jn4.i  ,  , 


h  -  W.s 

I  r^.  .  -  ^  fc  1 


Thus  the  new  general  filter  is  obtained  by  setting 

BO 

(  ^K4 1 , 0 , '  ’ ' ,  "I'nM  ,»i , '|n+ >  ,rH-t  ] 

*  [  f n,c  r  "  I  t«.‘t .  0 1  jn (  ^n+' ,n^ /  " ,  i ,  b o  ] 

Thus  the  general  filter  polynomial  satisfies  the  recursion 


In  particular:  |,(2>  =  l-.CZ)  + b.li) 

|ia)={,(Z>+-fe}.2^b,4) 

f,C2^=fW  +  -fe^2Hjd) 

BO  ^ 

{,(2)  •  ^.tz)  b-<2)  +  *ii  2^^  +  6^i2^t>3‘2  ) 


and  hence: 


I^CZ)  =  f.tz1 +4^,2  b,(^'/ +  t^,2*b,{j) +  -'l-^|,n-|Z''b,(i\ 


where  boC2)  =  bo  , 

are  constant  rntrice 


Q^cz)^ac  ,  b5(^)'bo 


SEI3MIC  ARRAYii 


A  seismic  disturbance  except  In  special  cases  does 
not  reacn  each  detector  in  an  array  at  the  same  instant^  and 
so  there  Is  a  time  difference  between  corresponding  phase  points 
of  the  disturbance  recorded  at  different  spatial  locations. 

This  time  difference  can  be  expr^^ssed  in  units  of  time  per  unit 
of  distance,  and  thus  is  the  reciprocal  of  the  plmse  velocity. 


Let  us  consider  the  following  model.  The  signals 
from  distant  events  arrive  at  the  array  with  horizontal  velo¬ 
cities  ranging  from  i.nfinity  (that  is,  the  case  of  vertical 
Incidence)  to  some  finite  minimum  value.  The  noise,  on  the 
other  hand,  is  assumed  to  have  horizontal  velocities  smaller 
than  the  minimum  set  for  the  signals,  and  such  lower  veloci¬ 
ties  are  indeed  to  be  expected  of  surface  wave  noise.  Thus 
in  this  model  seismic  si^ials  and  seismic  noise  have  different 
velocity  distributions,  and  therefore  improvement  in  signal- 
to-noise  ratio  can  be  expected  through  the  use  cl  a  device  that 
discriminates  against  noise  velocities.  Such  a  device  is 
called  a  velocity  filter,  and  for  the  model  that  we  have 
assumed  it  should  be  a  high-pass  velocity  filter,  that  is,  one 
that  passes  high  velocities  and  stops  low  velocities. 


Let  us  consider  a  plane  wave  front.  Because  a  two- 
dimensional  group  of  detectors  is  equivalent  to  a  one- 
dimensional  group  in  a  plane  perpendicular  to  the  wave  front. 

We  shall  limit  the  present  discussion  to  one  dimensional  arrays . 
A  plane  wave  of  a  fixed  temporal  frequency  ft  and  spatial 
frequency  jx  nay  written 


S  ww.  (.  2-nf,t+nTjx^  +  b  j 


all 


•  r  sn  r~.  ^  __  -  __  .  .  .  , 

-wii^vaiiL  Ff^OX’^Sf^rir  1  r-.rF 


u  X  5  y  xti  c  n  t- 


ig  a  pna 


pi  ia  o  e 


^  / 


/ft 


iS  the  perlud  f 


ci  F  C 


the 


^4  4?;  rXj_  4  s3  V 


CT  X  i  f  jii  r%  V'.  T  r 


5  ’^5  #=®  nP  •_»  O  t  r  »--  ^  ^  s  ^  .  -» 

vs^Aocity  ail 


v- 


.il- 

lx 


We  my  suppose  that  the  abov 
Units  as  rollows: 


fc  -jUantlties  are  expressed  In 


I't  =*  temporal  f  reQuenc''^  i  •-’  ^  i 

*Ae-iuenoj,,  1^4  cvcles  per  second  (c 

*/|t  =  period,  in  seconds  per  cycle  (s/c) 

=  spatial  frequency,  in  cycles  per  kilometer 

(c/kon) 

*/|x  =  Wavelength,  In  kilometers  per  cycle  (km/c) 

V  =  velocity.  In  kilometers  per  second  (km/s) 


■fyp  -d/jj/d/^tl  -  (wavelength)  /  (period) 


of  h  +  i  r*  ^  so-dalled  uniformly  distributed  array 

outbut,  ,  ^  to  the  sum  of  the 

outputs  from  _  _  -  -  ^ 

vi^i^  .  oec  l  igureD* 


Figure  b.  Uniformly  distributed  array  fh.7) 


Henct;  the  uutput  to  the  above  sinusoidal  input  wou; 
(normalized  by  the  factor  j/r\4i  ) 


1141  ^-0  ft  < 


which  is  equal  to 


im(  ^  5;.(2n^tt+TTb^+p) 


Cn.+o 

'  yi  ^ 


Thus  the  transfer  function,  namely  the  ratio  of  output  to  input, 
is  at  the  center  of  the  array 

If  we  let  ,  then  Ah.  becomes  the  transfer  ratio  of  a 

continuous  detector  whose  length  is  A  ;  we  have 

A  _  Sv>^TI-^x^ 

Aoo  “  TTjxA 

Let  us  Introduce  the  dimensionless  spatial  frequency  variable 
•  Then  we  have 


^;n(^nv) 

tn+()4!n(^) 


<u^^n/ 

irp 


The  figure  shows  graphs  of 


/  j(V)  ,  and  /\^iV)  . 


Figure  6.  Graphs  of  spatial  frequency  transfer  functions 


Provided  the  noise  corresponding  to  values  of  V  >6  is  neglige 
ible.  It  is  seen  that  increasing  the  number. ri-tt  ,  of  seismo¬ 
meters  in  the  array  beyond  10  can  yield  no  appreciable  improve¬ 
ment  in  velocity  discrimination  under  the  scheme  of  uniform 
spacing  and  uniform  weighting  of  the  seismometers  in  the  array. 

For  example,  let  us  suppose  that  the  ambient  noise  is 
spatially  organized  and  is  propagating  with  an  apparent  horizontal 
velocity  of  4  km/s .  Suppose  that  one  component  of  the  noise  has 
a  temporal  frequency  -ft  =  1  c /%]  and  that  the  other  component 
of  the  noise  has  a  temporal  frequency  ft  =*  0.5  c/s.  Suppose 
that  one  component  of  the  signal  Is  made  up  of  mantle  P-waves 
travelling  at  a  velocity  of  8  km/s  and  with  a  tempoial  frequency 
1  c/s;  that  the  other  component  of  the  signal  has  the  same 
velocity  (8  km/s)  and  a  temporal  frequency  ft  “  2  c/s.  Thus  we 
have  the  table  (where  wavelength  X  -  the  reciprocal  of 

the  spatial  frequency). 


TEMPOfiAL  / 

FREQUENCY 

NOISE  1 

Velocity  1 

y-  *  4  km/s  j 

Q  -  O' 5  % 

k-  I'O 

X  =  &  +.W ,  *  Y 

X=  f  iffi,  1^=  1 

.  —  -  ..  -  _  - -  ■  —  ■ 

i 

i 

! 

X  =  8iM,  1/*=^  i 

X  =  ^ 

TABI£:  Wavelength  X  and  dimer.s  ion  less  spatial  frequency  2/ 
vs.  velocity  V  and  temporal  frequency  ft  ,  where 
t  ®hd  A  -  4  km. 


Figure  7.  Transfer  function  of  AwCVJ, 

We  see  that  the  maximum  wavelength  that  can  be  completely 
filtered  out  is  equal  to  the  array  of  length  A  ,  so  the 
minimum  length  of  the  array  is  dictated  by  the  n^ximum  wave 
length  which  is  desirable  to  have  filtered  out.  Thus  we  have 
the  filtering  action  of  the  array  given  by  the  table: 


SIGNAL  1  NOISE  j 

V=  8  1  y,  4 

i 

i 

i 

H — - — - — i 

it  ■  ^*5  c/s  i 

(PASSED)  1 

j 

i  it  -  i  c/s 
i  (PASSED) 

- -  1  ■  '-j 

j 

it  -  1  c/s  j 

(STOPPED)  1 

f ^  2  c/S 

(STOPPED) 

j 

- - — - - - — _ _ 1 

TABLE:  PASS  and  STOP  bands  of  array,  with 

spread  A  »  ^  km. 
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Thus  for  8  kra/s  P-waves  this  array  with  uniform  summation  pro¬ 
vides  maximum  effectiveness  only  over  a  small  fraction  of  an 
octave  around  the  temporal  frequency  -Jt  =  1  c/s. 

Now  le  ;  us  consider  an  array  with  equally  spaced  de¬ 
tectors  but  where  each  detector  is  wel'5hted  by  an  arbitrary 
factor  before  sunsnatlon.  We  suppose  that  we  have  2h+l  detec¬ 
tors  of  arbitrary  sensitivity  at  locations  0,  ’ ,+A 

with  sensitivities  I,  Cli,  (Xi,"*  , O13  .  .  If  the  input 

is 

Cos 

then  the  output  is 

Co5t2Ttftt^f)(l+2i  CljC.s(2Tric-^)j 

so  the  transfer  ratio  is 

A  =  1+2^ 

flu  3=1  ^ 

We  see  that  is  a  periodic  Junction  of  .  When 

we  have 

t  ^ 

when  Tx  "Tf  we  have 

k  I  1 

Thus  the  period  of  Mk  corresponds  to  *  which  is  a 

wavelength  t  namely  the  detector  spacing  of  the  array. 

In  other  words,  the  array  cannot  distinguish  between  the 
following  two  situations: 
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Figure  8.  Illustration  of  aliasing. 


This  phenomenon  is  well-known,  ana  is  called  aliasing.  Thus  it 
is  desirable  to  choose  the  detector  spacing  such  that  the  100^ 
response  peaks  that  occur  after  ^>^*0  do  not  correspond  to 
wavelengths  that  are  associated  with  substantial  energy. 


Brustad, 


Suppose  the  following  values  of  (Xj  are  chosen 


and  Sider,  1958): 

Oo  =  1.000 

di  »  0,987 

at  =  0.592 

Oj  -  0.947 

a-j  =  0.473 

Oij  =.  0.888 

Qj  =  0.348 

cu  -  0.809 

0.224 

=  0.704 

0.105 

(Savit, 


The  transfer  ratio  in  terms  of  dimensionless  spatial  frequency 

j/  =  jjcA 

Graphs  of  (1^  and  Ak(V)  are  shown  In  the  figure 


Figure  9.  Weighting  function  Ot-T  and 
function  ^ 


corresponding  transfer 
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5.  VELOCITY  PILTERIjgC 

The  concept  of  filtering  seismic  signals  based  on 
their  apparent  velocities  is  treated  by  Texas  Instruments 
Staff  (1961),  Embree  et  al  (I963)  and  Fall  and  Graw  (I963). 
Such  velocity  filtering  for  the  nuclear  surveillance  problem 
is  implemented  by  j.ie  use  of  large  arrays  of  seismometers. 
Velocity  filtering,  then,  is  a  multichannel  filtering  process. 
The  process  uses  a  number  of  input  channels  from  the  seismo¬ 
meters  to  produce  each  output,  but  differs  from  conventional 
mixing  and  filtering  techniques  in  that  each  input  at  a 
different  delay  is  filtered  through  a  different  appropriately 
designed  filter  response  before  being  summed  to  form  the 
output. 


In  order  to  design  such  '  multichannel  filter,  a 
working  mo(?>l  of  signal  and  noise  must  be  formulated  in  the 
temporal  and  spatial  frequency  domain.  In  addition,  criteria 
for  Judging  the  performance  of  the  arrays  must  be  established. 
The  theoretical  requirements  for  optimum  processing  liave  to  be 
developed  with  the  quantitative  evaluation  of  the  performance. 

The  primary  objective  of  large  array  systems  is  to 
provide  a  better  picture  of  the  seismic  signals  when  they  are 
masked  by  ambient  noise.  A  large  array  should  make  it  possible 
to  detect  and  Identify  seismic  events  that  are  obscured  by  am¬ 
bient  noise  when  only  a  single  sensor  is  used. 

The  complexity  of  the  seismic  waves  that  arrive  at 
a  detection  station  is  well-known.  For  either  an  earthquake 
or  an  explosion  source,  the  following  sc^quence  of  seismic 
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wavej  is  produced: 


(^) 


longitudinal  body  waves . 
abon-^r?  velocity  of 

medium!  ^  depending  on  the 

transverse  body  waves, 
of^p-wlvll!  velocity 

horizontal  transverse 
surface  waves.  They  travel  along 

4  ^  velocity  of 

Rayleigh  waves,  or  vertical  surface 
waves  that  perform  a  retrograde 
elliptical  motion  with  its  plane 
lying  along  the  axis  of  propagation. 

The  P  and  s  waves  are  subject  to  reflections  and  refractions, 
and  the  surface  waves  undergo  dispersion.  The  result  is  a 
complex  sequence  of  oscillations  whose  nature  Is  governed  by 
the  Characteristics  of  the  transmission  paths  and  very  little 
by  the  source  Itself,  since  the  only  part  of  the  wave  train 
that  arrives  undisturbed  by  later  modes  Is  the  P  wave  (the 
first  arrival),  it  is  the  most  Important  signal  currently 
used  for  location  and  Identification  of  the  source.  Secondary 

Objectives  are  the  extraction  of  s-waves  and  surface  wa.es  In 
the  presence  of  noise. 

The  noise  obscuring  the  desired  signals  may  be 
categorized  as  follows; 


(a)  Ambient  coherent  microseismic  noise 
propagating  primarily  as  Rayleieh  ^ 
waves,  with  an  apparent  horizontal 
velocity  of  ?.5  to  3.5  km/s .  The 
isotropic  assumption  means  that  the 
noise  is  equally  likely  to  come  from 
any  azimuth. 
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(b)  Local  coherent  noise  from  sources 
such  as  factories,  lakes,  railroads, 
highways,  propagating  primarily  as 
Rayleigh  waves,  but  with  a  particular 
azimuth  of  propagation. 

(c)  Incoherent  noise  from  sources  or 
scatterers  within  the  array  area, 
including  locally  generated  wind 
noise , 

(d)  Incoherent  instrument  noise. 

The  P-wave  signal  may  be  characterized  as  follows: 

ia)  Coherent 

bJ  Equally  likely  to  come  from  any  azimuth 
c)  Apparent  horizontal  propagation  velocity 
of  8  to  15  km/s. 


The  signal  and  noise  can  be  shown  in  the  three 
dimensional  temporal  and  spatial  frequency  domain,  whose  axis 
are  example,  a  propagating  sinusoidal  plane 

wave  is  a  point  in  the  jt,  domain,  even  as  a  sine 

wave  6‘m  iTT-Jt't  is  a  point  in  the  -Jt  domain.  The  signal 
and  coherent  noise  in  the  working  model  are  continued  in  the 
conical  boundaries  shown  in  the  figure: 


Figure  10,  Signal  and  noise  in  three  dimensional  frequency  domain. 
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Assuming  that  the  signal  propa^tes  with  a  speed  of  from  8  to 
15  km/s,  and  is  equally-likely  to  arrive  from  any  azimuth,  the 
signal  power  lies  within  the  conical  boundaries  shown  in  the 
figure.  The  total  power  at  any  frequency  be  obtained 

by  integrating  over  all  'i'i-  for  that  -^t  >  and  It  is 
assumed  to  be  the  curve  shown  in  the  figure: 


/  X  ^  ^ 

7^  VH  C/A» 

Figure  11.  Signal  power  density  spectrum 

The  noise  power  lies  between  the  two  cones  corresponding  to 
velocities  of  2.5  to  3.5  km/s,  and  by  integrating  over 
and  assume  the  power  density  has  the  form  in  the  -ft 

domain  as  shown  in  the  figure:. 


ft  tff  C/U, 


Figure  12,  Coherent  noise  power  density  spectrum 
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Incoherent  noise  at  a  particular  frequency  is  uniformly 

distributed  over  the  entire  plane. 

For  the  Isotropic  model,  a  vertical  cut  in  the 
-^t>  \y  space  is  independent  of  azimuth,  as  shown  in 
the  figure: 


Figure  13 •  A  vertical  slice  of  Figure  10 


For  a  directional  signal  model,  the  three  dimensional  space 
must  be  retained. 


This  working  model  of  seismic  signal  and  noise  has 
the  following  important  features: 

(a)  The  signal  and  coherent  noise  are  com¬ 
pletely  separable  in  ft;  space. 

This  separability  is  unique  to  the 
seismic  case;  3n  radar,  sonar,  and  radio 
astronomy  the  signal  and  noise  propa¬ 
gation  velocities  are  identical  and  so 
separation  can  be  achieved  only  on  a 
directional  basis.  The  maximum,  improve¬ 
ment  in  such  cases  is  /n  (in  amplitude) 
in  signal-to-isotropic  noise  ratio  by  an 
H  element  array;  because  the  signal 
and  coherent  noise  are  completely  sepa¬ 
rable  in  the  seismic  case,  greater 
gains  are  possible,  and  the  achievable 
gain  is  dependent  on  the  temporal  and 
spatial  distribution  of  the  noise. 


co.pXe\®2!y 

Sfin  a??a"r  «ffectl5rij^re1iced 


1  <a.  ^  ^  consider  the  actual  implementation  of  the 

eloolty  filter.  We  want  the  filter  to  paaa  waves  In  the  vel- 

oolty  range  from  -V  to  +V  ,  so  the  desired  response  In  tem- 
poral  and  spatial  frequency  is 

I  0 

Afft,ix)=\  I  -jftv  ^  lit! 

V  Y 

0  -for  *^<|x 

This  transfer  function  is  depicted  in  the  diagram; 


Plgiire  14,  Desired  spectral 


response  function  of  velocity  filter 


At  a  given  frequency  jt  the  spatial  impulse 


response  is 


4 

Tx  Y 


/ 


dt. 


iM 


/.iil 

V 

I 

ITTX 


CosHTT^xX 


ik 


5tK  2Tr 


LW 

V 


TDi 


If  lftl“V  then  a(it,)i;)=~ 


5m21TX 


,  which  looks  like: 
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whereas  if  Htl=iV  then  ,  which  looks 

like: 


Figure  16.  Spatial  impulse  response  of  desired  filter  for 

11,1= iv; 

Therefore,  in  terms  of  arrays,  such  a  process  appears  at  low 
frequencies  as  a  very  long  array,  which  passes  long  wavelengths 
and  stops  short  wavelengths,  or  in  other  words  it  acts  as  a 
narrow  band  low  pass  filter  in  spatial  frequencies.  Such  a 
prcv'  ;88  appears  at  high  frequencies  as  a  very  short  array, 
which  passes  shorter  wavelengths,  or  in  other  words,  it  acts 
as  a  wider  band  low  pass  filter  in  spatial  frequencies. 


The  impulse  response  in  time  t  and  space  X  is 


given  by 


a(t 


4  .IM 

’  --  L-f 

=  (X(j^,x)  dijt 
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In  the  usual  case,  the  array  will  be  made  up  of  equally  spaced 
sensors  with  spacing  AX  ,  as  d*'  pic  ted  by: 


U— Xm 

J _ A _ _ A 


3^ 


I 


Figure  17*  Sensor  spacing 


We  suppose  tiiat  there  is  an  even  number  of  sensors,  and  we 
let  Xm,  denot(  the  distance  of  the  sensor  from  the 

center  of  the  array.  Thus  =  etc. 

The  process  may  be  implemented  by  passing  the  output  of  each 
sensor  through  a  temporal  frequency  filter  with  transfer 
function  specified  by 


<  n-rriM.  Y 


and  then  summing  the  outputs  of  all  the  temporal  frequency 
filters.  The  frequency  spectrum  CXC-^t.Xw)  looks  liic 


Figure  18,  Frequency  spectrum  of  (Xl-JtiXm) 
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At  any  given  temporal  frequency,  no  array  with  equally  spaced 
sensors  can  distinguish  between  zei‘o  cycles  per  2  ,  or 

any  other  integral  number  of  cycles  per  2iX  ,  as  shown  in 
the  figure: 


Figure  19.  Spatial  aliasing 


Here  =  -fx  is  the  Nyquist  spatial  frequency  (or  the 

folding  spatial  frequency),  and  the  equally ‘spaced  sensors 
result  in  a  periodic  response  in  spatial  frequency,  the  response 
being  repeated  at  the  interval 

“  AX 

and  the  desired  response  can  be  achieved  cmly  over  the  spatial 
frequency  band 
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At  temporal  frequencies  above 


V 

2  AX 


the  specified  pass  band  of  the  velocity  filter  extends  beyond 
the  Nyquist  spatial  frequency  »  so  the  temporal  frequency 

range  over  which  a  useful  response  may  be  obtained  is  limited 
to  temporal  frequencies  below  jx  .  Thus  a  sampling  Interval 
At  =  ■  y I  may  be  chosen,  and  no  information  is  lost  provided 

the  data  is  limited  to  frequencies  below  and  .If  we 

let  be  the  time  of  the  sample  point  in  the  is- 

crete  impulse  response  of  the  time-domain  filter,  then  we  ee 
that  V  satisfies 


— 

At 


The  required  coefficient  of  the 
domain  filter  for  the  sensor  is 


r 


time  point  of  the  time 


'n'Xm 
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1 


1l3[ 


j  Ajt  CwUirnAtjo  SUciTT^^At'jt) 


Now  the  following  integration  formula  occurs  in  standard 
tables: 


f  w  ax  Coj  bx  ix  = -^ 

j  ii  (X-b  ^  oi+b  J 

provided  .  Hence 


a(T..x«)=~c-4-) 


indt(^-n)  2iTAt(K* 


AX 


.4d) 


■{t=0 


Cirfir(-^-t-K)  j  I 

2Tt*t(.|a  w)  2Trati^-»)  2Tr«tl^+ii1 


and  since 


m-l 

hi-fX 


H  w>o 

m-<0  where  Wb  is  an  integer 
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and  since 


is  an  Integer, 

i 


it  follows  that 


is  equal  to  an  integer  divided  by  ^  ,  so  the  first  two  terms 

in  the  expression  for  acTn,X«)  above  vanish.  Hence 


(XCTn.Xw) 


1  ] 

ITT  At 

^  J 

I 


Xm 


-h 

AX  AX 


+(%-«) 


I 


I 


where 


Trt  =  KAt  5  n  =  0,d: 


'  (m-j)AX 

,  ori-i 


tn  =  i.2,3,'" 


It  is  seen  that  the  linear  operator  Q-CTiiXw)  is  symmetric 
about  the  origin  in  both  space  and  time. 


The  theoretical  frequency  response  can  be  recovered 


Atit.ix)  = 


2  I 

^=-00 

m;to 


by 


,.,,ll,niq|ii>||'l|il||li|||l<nnl!lllli:illllllllllllll|!ll|lll|l|lllll!'»"1lllln^  . . 
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In  practice  only  a  finite  nuntier  of  fl  and  Wl  values  can  be 
used,  so  the  theoretical  frequency  response  is  only  approximated. 
Incase  Tl*  0,±  and  1 1  ,±2,'-',±6  then  the 

actual  response  is  as  depicted  in  the  figure: 


fnqoetcy  AciutJ  {rt^ticy  ff^ponse 


Figure  20.  Theoretical  vs.  actual  frequency  response 


iliilllllliKll: 


riif  .  ^  illustrate  the  actual  process  of  velocity 

rigurl:”®' 


Hr* 


»5 


*  * 

ir 


1/^  ^ 
r- 


K 


v- 


<nmvr  or 

VUOCtTY 

AUfR 


I  '  t^-o./ 


JJ  Y  lOMm/jnu,  ^  cuf-aff 


^  f  t 

«65o!$* 

oi^HWs^  wwn  r^sses 


Figure  22.  Effeote  of  velocity  fUteri.ng  m  the  tl«e-spao-  domain 
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Let  us  choose  the  cut-off  velocity  V  equal  to  say  10  km/sec 
and  let  us  suppose  the  spacing  of  the  sensors  is  AX  =  1  km» 
Then  the  Nyquist  spatial  frequency  is 


1 

1  AX 


and  the  Nyquist  temporal  frequency  is 


AX 

At 


=  (i)(io)=  5 


Hence  the  discrete  time  spacing  is  chosen  to  be 


At  = 


_ I _ 

”1CS)”  10 


s<c  *=  o.l 


Hie  foregoing  figure  portrays  four  events  with  velocities  of 
5  km/s,  10  km/s  =  the  cut-off  velocity,  20  km/sec,  and  infinite, 
from  left  to  right.  The  waveform  is  a  sharp  transient  that  is 
identical  for  all  events.  Below  the  input  record  is  the  output 
of  the  velocity  filter  designed  to  pass  events  with  a  velocity 
greater  than  the  cut-off  velocity  of  10  km/sec.  The  output 
trace  corresponds  to  the  center  point  of  the  array.  Interpre¬ 
ting  the  output  trace,  we  observe  the  following: 


(a)  Events  within  the  pass  region,  that  is, 
with  velocity  greater  than  the  cut-off 
velocity  V  =•  10  km/sec  are  passed  by 
the  process  with  virtually  no  waveform 
distortion . 

(b)  Events  having  a  velocity  equal  to  the 
cut-off  velocity  V  »  10  km/sec  are 
attenuated  by  6  db*  as  shown  by  the 
actual  frequency  rtsponse  contour 
above,  but  still  with  virtually  no 
waveform  distortion. 

(c)  Events  having  a  velocity  less  than  the 
cut-off  velocity  are  generally  attenuated 
in  amplitude  by  20  db  or  more. 
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Velocity  filtering  offers  the  capability  to  enhance 
signal-to-noise  ratio  significantly  and  without  danger  of 
deteriorating  signal  waveform  as  no  signal  bandwidth  is 
sacrificed.  In  other  words,  it  makes  it  possible  to  process 
seismic  array  data  in  such  a  way  that  all  seismic  events  with 
velocities  in  a  given  range  are  preserved  with  no  alternation 
over  a  wide  frequency  band,  while  all  seismic  events  with 
velocities  outside  the  specified  range  are  uniformly  and 
severely  attenuated.  Velocity  filtering  may  be  looked  upon 
as  a  process  whereby  it  is  possible  to  combine  the  elements  of 
a  line  array  in  a  nanner  resulting  in  a  directed  beam  with  low 
side  lobes,  where  the  beam  width  and  side  lobes  are  essentially 
independent  of  frequency.  By  applying  the  velocity  filtering 
process  to  a  noisy  array,  an  output  may  be  obtained  that  has 
all  events  within  a  specifie'i  velocity  range  perfectly  pre¬ 
served  and  events  without  this  range  essentially  eliminated, 
a  result  which  is  impossible  by  conventional  array  usage. 

In  order  to  effectively  utilize  velocity  filtering 
in  practice,  a  noise  and  signal  analysis  procedure  imjst  be 
followed  that  will  furnish  the  needed  Information  for  effective 
array  design  and  emplacement.  Of  course,  it  is  feasible  to 
emplace  a  standard  array  layout  at  each  site  without  prior 
measurement  and  analysis  of  the  noise.  The  signal  and  noise 
parameters  for  the  emplaced  sensors  would  then  be  estimated 
from  the  actual  array  output,  and  then  these  i^raraeters  can  be 
used  in  setting  the  filter  responses  for  the  processor.  Thus 
if  a  standara  array  ?ayout  were  required  for  political  reasons, 
and  no  measurements  were  allowed  prior  to  site  selection,  it 
still  would  be  feasible  to  utilize  the  array  effectively. 
Nevertheless,  a  more  effective  network  can  be  established  by 
obtaining  local  noise  and  signal  characteristics  and  utilizing 
them  for  local  optimization  of  the  array  emplacements.  Tiius 
an  effective  noise  analysis  and  array  design  procedure  can  be 
worked  out  in  which  detailed  noise  measurements  are  made,  and 
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a  data  analysis  is  performed  to  establish  the  temporal  and 
spatial  frequency  characteristics.  The  necessary  array  length 
and  number  of  sensors  '*0  provide  good  slgnal-to-nolse  ratios  can 
then  be  determined  from  the  noise  characteristics.  The  per¬ 
formance  of  the  array  is  then  checked  out,  and  It  is  put  into 
production.  Nevertheless,  there  are  many  complications  to  such 
a  relatively  straight-forward  procedure. 

The  problem  of  noise  measurement  can  be  divided  into 
three  phases: 

(a)  Measurement  of  signal  and  noise  for  the 
purpose  of  selecting  the  site. 

(b)  Measurement  of  signal  and  noise  at  the 
selected  site  for  the  purpose  of  establish¬ 
ing  the  permanent  array  layout, 

(c)  Measurement  of  signal  and  noise  on  the 
emplaced  array  for  the  purpose  of  deter¬ 
mining  the  necessary  filter  settings. 


In  order  to  provide  the  desired  information  about 
tbf  noise  and  signal,  data  analyses  must  be  performed.  It  is 
anticipated  that  a  real-time  analysis  of  the  data  will  not  be 
possJbxei  thus  selected  samples  of  noise  and  signal  data  are 
the  inputs  to  a  digital  computer  analysis  program. 


The  parameters  of  the  most  Importance  in  the  design 
of  an  array  are  the  distribution  of  noise  and  signal  versus 
apparent  horizontal  propa^tion  velocity,  azimuth,  spatial 
frequency  spectrum,  temporal  frequency  spectrum,  time  stability 
of  the  noise  characteristics,  coherence  versus  frequency  and 
sensor  separation.  The  desired  information  can  be  derived 
from  the  data  by  operating  on  the  cross-correlation  function 
and  its  spectral  density.  Also  the  cross-power  spectra  of 
signal  and  noise  constitute  necessary  Information  for  the 
design  of  the  ideal  piocessor.  Thus  the  accuracy  in  the 
determination  of  the  necessary  parameters  is  dependent  upon 
the  accuracy  with  which  the  density  spectra  can  be  estimated, 
and  so  this  problem  will  now  be  discussed  in  some  detail. 


The  development  of  a  theory  of  spectral  analysis  of 
empirically  observed  signals  begins  with  the  introduction  of  the 
notion  of  the  sample  autocorrelation  function.  Suppose  we 
let 


denote  the  sampled  values  of  the  signal  of  interest, 
sample  autocorrelation  is 


(  -l(X,X,rt  +  XiXi+t  + - l-Xn-tXh) 


+,lt)  =  '  0 


Then  the 
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The  periodogram,  or  Fourier  transform  of  the 
correlation,  is 

iit‘h 


iample  auto- 


a 

2 

The  perlodogram  may  also  be  written  as 

=4’„i'»+2, 1  coswft  ■ 

I 


and  also 


=  i  CoaTTtf§„Ci^<l^. 


functltlr":*  autocorrelation.  Is  an  even 


ilm.  i 

K-»oO  J  ^ 


^ere  ^Xt)  is  the  autocorrelation  function  of  the  process 
Then  for  every  continuous  function  Gflp  we  have 


n^r^o  A  ^ = J!::,  i,  s-.f)  i 

*2 

~  j  5^)  Z 

=A 


^here 


7B 


is  the  spectral  density  of  the 
the  sample  autocorrelation  is 


process. 


The  expected  value  of 


so  its  expected  value  is 
the  true  autocorrelation 


biased  by  a  factor  of 


from 


autocorrelations  given  by  ™  ^  combination  of  the  true 

Now  the  expression  in  brackets  1h  -m  n 

assume  ttt «  Oi-t  of  the  w  <  n  ^'’ansf orm  (as  we 

CT  r  oi  the  Weights  QL>  •  « 

la  called  the  spect^l  window  A(f)  ),  ti;at  is 

OO 

XT***© 

I^t  US  now  assume  that  the  statinr,=  vv,r 
normal  with  zero  mean.  Then  s^tionary  process  is 


Act,)] 4  C0vr{  r 


where 


UiK) 

rt 


as  Yi  — ,  Thus 


as  Yi  - — >03 


But 


~2 

=  2  Cosairito 


Hence 


~ J  §\{)  e"’'*'^^'  2  c^saitn.)  ij 
2 


2  I  ^ 

~J  §  (4)  CoiC2njt,lCwairfts)dj. 


as  10,  — »oo  . 


Thus  the  deviations  of  the 
autocorrelation 


sample  autocorrelation  from  the  true 


is  approximately  normal  with  means 

E  { =  E  {  4cu}  -  4'C»  =  ( !  )  4>t^)  - 


0 


QC 


n 


8o 


and  covariances 


C5V I  Snitii  j  “  Co^  ^  i 


2. 

^j  §\()Cssnijt,  cwaitjtitilj- 


The  spectral  band  power  estimate  is  an  estimate  of 
tho  spectral  mass  in  the  interval  which  is 

I  i(jidi+J  §c4)dj] 


The  spectral  window  is 


-Jor  -  ^  ^  ^ 


which  is  an  even  function  of  X  rnK^ 

oi  y  ,  The  corresponding  weights 

are  1 

Qt=f  =2^  Ai|)  CoS  nr^t  df 


3Tf  C^.HH  -  jin  I 

iTTt  *lr 


to^lJrlot  5JH2ir-R\ 
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Thus,  in  order  to  estimate 


iat4rt)=  I  ■  j 

-1 


we  could  use  the  estimate  (the  spectral  band  power  estimate) 


jL 


and  since  -0  for  jvl^K,  this  estimate  can  be 

written 

I  ati't) »  I  'H  =  j 


ts'OO 


In  other  words,  the  spectral  band  power  estimate  is  equal  to 
the  energy  in  the  periodogram  in  the  spectral  band  of  interest. 
The  expected  value  of  this  estimate  is 


=  I  Qt 


/L  V 
ltl<K 


=  I  0tWi4’wih 


<--*0 

Now  the  function  defined  by 


-{or  ftl^n 
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serves  as  a  weighting  function  to  improve  the  convergence  of 

I  at 

to  as  V\“-^Qo  .  In  other  words 

1  QtUtCtsMt 

HKH 

is  a  better  approximation  to  AC^)  .  The  weighting  function 
UJt  Just  given  is  the  well-known  Pejer  weighting  function; 
of  course,  any  number  of  other  choices  are  possible. 


by 


Then 


For  example,  suppose  we  estimate  the  autocorrelation 


0  'foir  }r>Kt 


^ov  1*^1  >H| 


so  the  expected  value  of  the  spectral  band  estimate  is 

e(  I  =  I  atilt)  - 1  QtiJtitt) 


Where 


/ 


I 


I  0 


4®*“  Itj  ^ 
hr  |t|>n,. 


This  is  the  well-known  Dirichlet  weighting  function. 
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The  statistic  ^  under  the  stated 

conditions  is  normal  with  mean 

1 

\m  k 


and  variance 

lt,Kn 

ltj<n  jbl^n  i/l 


Thus  the  above  statistic  suffers  from  two  kinds  of  errors.  One, 
the  bias  comes  up  because  the  mean  of  the  statistic  is  not  equal 
to 

\ 

t=-!»  4 


8k 


but  is  equal  to 


|t(<a 


f  1  Qt  i^rt  Col  2Hf  §(  dl  j 
l«<n 


^rthor  dlscuaslon  of  these  points  is  given  In  the  work  of 
^rzen  (I96I,  1964)  and  Blackman  and  Tukey  (I958). 
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7.  DISCUSSION  OF  "PECTRAL  WINDOWS 


As  we  have  seen,  narrow  time  signals  have  spectra 
that  are  wide,  and  conversely.  In  order  to  give  a  quantitative 
measure  to  this  observation,  we  need  to  define  the  durations 
of  a  signal  and  its  Fourier  transform  in  a  simple  and  useful 
way,  but  no  single  definition  can  be  suitable  for  all  possible 
signals . 


For  a  measure  of  the  duration  of  a  time  signal  the  second 

moment  of  about  its  mean  or  some  other  suitably  chosen 

point  can  be  used,  which  we  can  take  to  be  the  origin  t»0 
From  we  obtain  the  Fourier  transform  pair 


where 


~  .  By  Parseval's  formula  we  have 


We  therefore  conclude  that  high  ripple  in  the  amplitude  A  or 
phase  4*  results  in  signals  of  long  time  duration.  Among  all 
functions  with  the  same  amplitude  minimum  phase  one 

has  the  shortest  time  duration. 
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To  simplify  notation,  let  us  assume  that  the  energy 
of  the  signals  under  consideration  equals  one; 


The  Schwarz  inequality  is 


where  the  two  sides  are  equal  only  if 
to  .  Let  us  insert  into  this  inequality 


is  proportional 
the  functions 


Hence 


Integrating  the  left  hand  side  Integral  by  parts  we  have 

^QO 

where  we  have  assumed  I  is  real  and  !§  - 
Prom  it  follows  that 


0  as  t— >iOo 
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Hence 


VfAi 

X  J  ^LoC 

If  we  let  ~  J .'  and  then 


we  have 

The  equality  sign  holds  if  where 

Hk  is  a  constant.  Solving,  we  obtain 

=  C 

Choosing  "fe  ■=  ■“  2  ai  “  “ ^  we  obtain 


^ct:i==C 

The  transform  of  is 


Git)-CJ3p@e-i>f''"*''.fipc;5ne 


_i3ik* 


Thus  a  narrow  pulse  has  a  wide  spectrum  Q{^)  ,  and 

conversely.  This  pulse  ^it)  is  the  well-known  Gaussian 


Dulse. 


Another  set  of  well-known  pulses  are; 


j 


for  various  values  of  Yi,  .  These  are  called  the  Ittit  power 

cosine  pulses.  The  transform  of  is 

't 


G(b=  f  Jt. 

y'tr 
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Let  us  choose  our  units  so  "^=2  ,  so 

|or  |t|  4  j 

0  |i|>|. 

Then 

Qc{)*=  /  Coi^H  dt 


I 

CoS  iTift  dt 


A  HJ  CoiWf 

^-r~"  - -  wH«»i  n=i,3.5,T,-' 

|l»o 


=  2  _ ^’'n 

-k*l 


volvih  n=-i,4.^ 
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For  example,  if 

and 

If  on  the  other 

and 

Gi  = 

Letting  it 

Gc4>  = 


Tl  *  I  then 
toiJ^ 

4  tollTTf 

IT  1-Ib4^ 

hand  we  had  let  t  -  I  (and  yv.=  j  )  then 

'  tojiTt 

.  0  W>i 

1 
/  2 

C.oiTrt  CoS^^^t  dt 

we  have 

to5f  C«^2TTj|  dt' 


(t|  ^  ! 

(  "t  -i) 


1  4 

2  ^ 


from  the  foregoing  reaiilt. 


90 


Hence 


r  lis  -  A  Coi^i 


(  't-\) 


Now  suppose  yi  =•  2  .  Then 


CoS^-^  =  -^  +  ^  tos 


iTTt 


0 


Itl^l 

K  ltl>|. 


This  is  the  well-known  pulse  named  after  the  Austrian  meteor¬ 
ologist  Julius  von  Hann,  and  called  “hanning"  by  Blackman  and 
Tukey  (1958). 

If  we  let  ,  then  the  von  Hann  pulse  becomes 

^lt)  =  l(l+Co>^)  -ior  Itl^T 


Now  may  be  written 


1  I  I 

5ct)  =  j(i  +  ie  ^+ie  T  ) 


Httct) 


where  ,  read  “rectangle  of  t  “  is  defined  as 


(  \ 


(  G 


|t|4T 
-Jor  |t|  >T. 


The  transform  of  r^ctC't)  ig 


=  RECT(f)  *=  )  Ct5llT4t  it 

-T 


-2T 


i'rriT 


is 


Hence  the  transform  of 

\  Slnl^T  .  \  6:ftlTr(f^  ^+±  lirCf  jr^"^ 
'‘T  TTt  ■’■? 


imlTV^T  f 

t  t  i 

_ 

“  1  1 

TTi  iir(H^>  :i-Tr(i  J 

_  5in  iMT 

1 

1 

- 

2 

1  -i 

=  i’.nlTtjT 


-Jt 

_ r _ 

nij(caTr})*-(.i)h 


_  /■jfi-T.')  _ ! _ 

2TtfT  ^  TT’-cnrlT)^ 


If  we  let  0l  =  ITTiT  ,  then 


Gi4)=TT'T 


5tn-Qt 

CH 


I 

TT^-d^  ‘ 


■H'littniiiimtuitti 
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The  spectrum  is  thus  the 

fn  a  ^  function,  correspondims 

the  “ 

.  When  (i=±lr  ,  this  factor  is 

aTT.%" e  also  vanlohes 

ac  a»xT}  .  'j’j^e  first 

n-  ,  51  thus  suppressed 

Fo^easing  values  Of  N,  .  greater  than  TT  ,  the  factor 

decreases  rapidly,  thus  reducing  the  Intensity 

lobes.  This  Is  of  course  why  a  cosine-squared 
pulse  is  used  in  Dreferenfa  ^-k  a. 

estiman^.  rT  f  ‘^he  rectangular  pulse  In  spectral 

estimation.  The  following  figure  shows  Qglj)  for  the  rect» 

lar  pulse  for  jtUT  and  =0  for  |t|  >T  and 

the  cosine-squared  pulse;  ^ 


Figure  23.  spectra  of  rectangular  and  cosme-squared  pulses 
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8,  REUVnONSHiP  BETWEEN  CONTINUOUS  Afffi  DISCRETE  APERTURE 

fun:tions 


Suppose  that  we  have  an  array  of  equidistant  elements 
C  AX-b  *=  d)  connected  in  phase  with  weighting  T (X^  )  , 

See  figure  24. 


1 

1 

i 

i 

1 

1 _ 

1 

T 

3 

\fd-n 

Figxore  24.  A  six-element  array 


We  shall  only  discuss  the  symmetrical  array  consisting  of  2ix 
elements!  a  similar  theory  can  be  developed  for  an  odd  number 
of  elements , 


The  pattern  of  this  type  of  array  as  a  function  of 
direction  &  is  given  by 


llr  Xi^ 
X 


Let  us  indroduce  the  new  variables 


M  _  t'fc’iidi  _ft  I  >  I 


ihdL  2n-dl 
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and 


Ind 


Then  we  have 

|A.KC^.h2|XT*c..2Tr)c»5 

The  function  |  is  periodic  with  period  2h.  , 

Let  U3  now  consider  a  contin^ ‘^us  aperture  function 
Tc^^)  of  width  L*^Hd  passing  through  the  points  (X'jt.Tfc) 
of  the  array.  The  pattern  due  to  this  continuous  function  is 

A^f5)=  j"*  Tc(>:)Co5  2Trx5  dx 

CX 

2 

in  comparison  to  the  pattern  of  the  array,  which  we  recall  was 

In  order  to  show  the  relation  between  and  Ac(5) 

we  represent  the  array  T^,  as  the  product  of  with  a 

uniform  array  of  infinite  extent  and  reduced  Interval  1/1k 

J.  =  eombj^  * 

^  In. 


Here  we  use  the  notation 


CombpUcx)=  I  UCi^) 

See  Figure  25: 


2/n 


Figure  25 ^  T*  as  the  product  of  Tc^)  with  a  comb  function 


As  a  result,  the  spectrum  of  T*  is  the  convolution  of  the 
spectra  of  tombX  and  Tc<‘^)  *  Because  the  spectrum  of 

Comb-^  Is  another  uniform  array  of  Infinite  extent  with 

period  2itt«  t/di)  we  have 

Ae5' 


or 
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where  is  defined  by 

Oo 

ThU  is  the  well-known  result  of  aliasing  (Blackman  and  Tukey, 
iypo#  Jacquinot  and  Roizen-Dossler,  1964), 


9. 


LISTIMJ  OP  VARIOUS  APERTURE  FUNCTIONS  AND  THEIR  WINDOWS 


The  form  of  an  aperture  function  is  indeterminate  to 
the  extent  of  a  proportionality  factor.  In  order  to  compare 
different  functions  some  convention  must  be  chosen  to  fix 
this  factor.  We  suppose  that  the  real  aperture  function 
vanishes  outside  the  interval  ,  We  want  the 

function  to  have  the  maximum  possible  transmission.  Because  the 
Intensity  at  the  center  of  the  pattern  produced  by  the  uniform 
aperture  function  for  is 

=  [G.toi]  “  [  j,  ]  =  I  . 

‘5 


we  shall  consider  in  the  case  of  a  non-unifurm  aperture  function 
the  ratio 


1(0) 

1.(0) 


We  wish  to  decrease  side-lobes  as  much  as  possible  without 
decreasing  1^0)>/l,(0)  too  much. 


The  total  energy  transmitted  by  the  aperture  is 
proportional  to 

■t'  = 

-t 

which  for  a  uniform  aperture  function  is 

Therefore  we  are  concerned  by  the  ratio 
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Our  object  must  be  to  decre'^se  the  side-lobe  level  without 
decreasing  too  much. 

The  width  of  the  main  lobe  of  the  pattern  is 

also  of  interest.  Let  us  as  a  convention  put  the  width  A|o 
of  the  nain  lobe  of  the  pattern  from  the  uniform  aperture 
function  equal  to  unity;  this  unit  width  turns  out  to  be 
the  distance  between  the  two  points  with  ordinate  0.4o5.  We 
then  extend  this  definition  to  other  patterns;  their  widths  Al 
will  be  taken  always  to  be  the  distance  between  points  of  the 
normalized  ordinate  0,405. 


Aperture 


To  confirm  our  notation  in  this  section, 

0 


Pattern  ; 
(or  window) 


we  have: 

!>f| 


Intensity 


1(6)  *  i 


Zero  subscripts  "efer  to  the  uniform  aperture  function  3 

for  .  It  turns  out  also  that  all  the  aperture  functions 

we  consider  are  even  functions  of  X  so  <»*>  is  real. 
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1.  Triangular  aperture:  ^iX)=:l“j2x 


1(8^)  _  $in*^  1  “^8 

I'O)  “  (Yite)**- 


M 

I.W) 


=  0.15  , 


At 

Al!« 


l.4fe 


2.  Sinus  cardinal  aperture; 


1(6)  . 
1(0)  ' 


IfiTT  j 


1.- 

t. 


6iin 

TT 


=  0.45, 


I^lC)  V  IT  ) 


0,15  y 


aH 

^ic 


1.42 


3.  Cosine  apervure:  ^tl)  =  CcilTX 

1(6)  _  Co>^Tr6 

1(0)  “  (i-H-6")^ 


M=iL=o4o,  4. 

1,(0)  Tr‘  At 


10-0 


^  I  / 

4 ,  Gaussian  aperture:  ^OC)  ~  ^  ^ 


(Note:  This  It®)/lt0)  does  not  represent  the  intensity 
exactly,  as  it  neglects  the  cut-off  of  ^tx)  at  *) 


3l 


5.  (l-'fX** 

aperture: 

^tx)  =  U-4x^)^ 

ne)  _  1 

8  JiOTO)  'I 

2 

l<.o)  1 

istO.H-S 

) 

MjlILf.0,31 

lAB 

t. 

Lio) 

'  Al. 

6.  {i-4xy 

aperture: 

^cx)  =  (l 

l-4x*f 

Iiei  _ ' 

KO 

[  A^t6>]  =  1 

[isjf 

1 

i 

2 

^=o.4t 

\ 

=  0.i8 . 

Ai  , 

=  1.54 

t; 

U6) 

/ 

Ai. 

Note: 


Here  wnere  Jy  is  iessel 

function  of  order  J/. 
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7.  (  \  ~  4’^  aperture: 


r(.V'^)i'\ 


-j^CT  V  ^ 
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